Nonuniversal scaling behavior of Barkhausen noise 



Bosiljka Tadic"^ 

Jozef Stefan Institute, University of Ljubljana, P.O. Box 3000, 1001-Ljubljana, Slovenia 



We simulate Barkhausen avalanches on fractal clusters in a two-dimensional diluted Ising ferro- 
magnet with an effective Gaussian random field. We vary the concentration of defect sites c and 
find a scaling region for moderate disorder, where the distribution of avalanche sizes has the form 
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The exponents r(c) for size and a(c) for length distribution, 



and the fractal dimension of avalanches -Ds(c) satisfy the scaling relation Ds{c)t{c) — q(c). For 
fixed disorder the exponents vary with driving rate in agreement with experiments on amorphous 
Si-Fe alloys. 



Barkhausen noise (BN) occurs at low temperatures 
when a disordered ferromagnetic sample is slowly driven 
by the external magnetic field. A small ramp in the field 
triggers one domain, and the perturbation spreads to 
the neighboring domains producing an avalanche, which 
manifests itself as a series of jumps in the magnetiza- 
tion as the system passes from one metastable state to 
another. Hysteresis in the field-magnetization plane is 
inherent to these processes. If the field is ramped up 
and down close to the coercive field of the sample, the 
avalanches show a broad distribution of sizes, although 
the triggering mechanisms are strictly local and fully de- 
terministic. The avalanche-like dynamics is controlled by 
the spatial distribution of quenched-in defects of various 
types, and is not influenced by the thermal fluctuations 
(as long as the system is below its ordering temperature) . 

Besides its practical applications , Barkhausen noise 
has been studied both experimentally 1^-^ and theo- 
retically as an interesting example of collective 
dynamical behavior far from equilibrium. Recent ex- 
perimental studies of BN were made on amorphous al- 
loys such as metglass 2605S-2 and alumel Q, metglass 
2606TCA 11, VITRO VAC 6025X @, and on Fe-Ni-Co 
and Si-Fe |^-^ alloys. Numerous data collected in 
these and other experiments show that the power spec- 
trum decays with frequency as ^ uj~'^ (1.5< cj) <2) and 
the distributions of size, duration, and energy associated 
with Barkhausen jumps exhibit a power-law behavior 
over a few decades with a cutoff. More detailed anal- 
ysis reveals that the measured critical exponents for 
various distributions obey certain scaling relations, which 
are also derived in Ref. jj] assuming a deflnite type of the 
elementary pulse. 

Experiments were typically done on one sample and no 
special care was taken in controlling the strength (and 
type) of disorder, which is responsible for avalanche-like 
dynamics. Other details which influence the statistics of 
BN were studied. In particular, it was found in Ref. ||] 
that in 1.8% Si-Fe alloy the exponents characterizing the 
distributions of size and duration of Barkhausen pulses 



appear to vary linearly with the driving rate of the ex- 
ternal field. Urbach et al. ^ studied the importance of 
demagnetizing effects. 

The dynamics of field-driven random fcrromagnets at 
zero temperature has been studied by numerical simula- 
tions using Ising model with random fields (RFIM) [p|jTc|] 
and with random bonds of the spin-glass type (RBIM) 
p^ , p^ . It has been recognized that in both cases a crit- 
ical level of disorder exists at which the response to the 
external field becomes slow, with a broad distribution 
of avalanche sizes. This disorder-driven phase transi- 
tion has been reported to have unusually large critical 
region — 50% above the critical point — in the case of 3- 
dimensional RFIM |o). 

In the present work we would like to point out the 
relevance of another type of defects to the scaling behav- 
ior of BN. In magnetic systems random fields are gener- 
ated on a semimacroscopic scale by other types of dis- 
order (i.e., random site or random bond defects) when 
the system is placed in an external magnetic field [Q. 
At the asymptotic thermally driven phase transition ran- 
dom fields dominate, leading to a new universality class 
of (equilibrium) phase transitions However, far from 
equilibrium where the Barkhausen noise is produced, the 
impurities of this type may still be relevant and are com- 
peting with the random fields in the pinning of domain 
boundaries. In particular, in the case of Si-Fe alloys the 
sample is diluted by nonmagnetic ions, whereas in the 
case of Fe-Ni-Co alloys and metglass samples the type of 
defects can be described by the random anisotropy and 
random-bond models, respectively. 

We consider the diluted ferromagnetic Ising model 
with concentration c of empty (i.e., zero spin) sites, and 
add a small Gaussian random field with variance / at 
each occupied site. It should be stressed that our model 
differs from the ones considered in Refs. ||9|-|l^ in two cru- 
cial aspects: (a) There is a larger parameter space in the 
(c, /)-plane, where the scaling behavior of the avalanches 
could emerge, and a well defined percolation threshold 
along the c-axis; (b) For c 7^ the avalanches contribut- 
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ing to Barkhausen noise develop on noncompact clusters 
of spins (see Fig. 1), the compactness of which is tuned 
via the parameter c. This leads to new prominent fea- 
tures of BN in diluted ferromagnets . 

We start with a (coarse-grained) Ising model 
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±1 and the interactions are Ji- 



tween nearst neighbor occupied sites, and J^, 



Jo be- 
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least one site is a defect. We use a 2-dimensional square 
lattice with a fraction c of randomly distributed defect 
(i.e., Si — 0) sites and periodic boundary conditions. As 
already mentioned, the local random fields hi are gener- 
ated on a semimacroscopic scale. We assume Gaussian 
distributed random fields with zero mean and variance 
/ at each occupied site. The disorder is considered as 
quenched. 

The system is driven by small ramps of the external 
field AH starting from a strong negative field —Hmax 
(for practical reasons Hmax < 4). Parallel updating of 
spin states is applied according to the rule Si[t + 1) = 
sgn /i'°^(<). For each rump of the external field the sys- 
tem is updated as long as there are no more unaligned 
spins. Clusters of flipped spins are monitored (an exam- 
ple is shown in Fig. 1 |17|) and the distribution of cluster 
sizes D{s) is determined. The distribution is integrated 
over the whole hysteresis loop and averaged over total 
number of samples, ranging up to 400. We first show the 
results for a fixed driving rate 8h = AH /Hmax = 0-01. 

In Fig. 2 the integrated probability distribution D{s) 
of avalanches of size s or larger is shown for fixed random- 
field variance / — 0.5 and for few values of concentration 
c (varied along the dashed line on the phase diagram 
in Fig. 3). For very weak disorder close to the zero- 
temperature ferromagnetic fixed point (lower right cor- 
ner in Fig. 3) the system is ordered and the external field 
fiips the magnetization via a first order phase transition. 
Consequently, the hysteresis loop is rectangular and the 
integrated probability distribution shows zero slope (top 
curve in Fig. 2). At fixed / = 0.5 a finite slope appears 
first for c « 0.05, increasing gradually with c in the range 
0.05<c<0.35. Assuming 
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and fitting the straight sections of the lines in Fig. 2 
to Eq. (^), we determine the exponent t(c) for 0.05 < 
c < 0.3. We find that the exponent 1 -I- t(c) varies 
from 1.155± 0.002 at c =0.05 to 1.373 ± 0.001 at c =0.3 
(see inset to Fig. 1). It should be noted that the loss 
of the power-law behavior at a finite cutoff size is not 
exponential, but rather a stretched exponential with a 
c-dependent exponent (see also Ref. Q). 

In order to study the cutoff systematically, we varied 
the size of the lattice L at several fixed values of c. For 



the integrated probability distribution D{s,c,L) the fol- 
lowing scaling form is appropriate | |T9[ | 
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where a is the scaling exponent for the probability distri- 
bution of avalanches of length £ or longer, and Dg is the 
fractal dimension of avalanches (not to be confused with 
the fractal dimension of percolation clusters). The expo- 
nent Ac is associated with the presence of defects. If the 
scaling function on the right hand side of Eq. (|^) depends 
explicitly on ci^" with Ac > 0, the criticality (and the as- 
sociated power-law behavior) could be maintained only 
if the disorder is kept at a critical value, say, cq, i.e., we 
are dealing with a (nonequilibrium) phase transition. For 
c 7^ Co the system is subcritical and the coherence length 
should be specified as a function of disorder ^(c — co) . On 
the other hand, the parameter c may appear implicitly 
via tuning of the exponents a and Dg in (^), while the 
simple finite-size scaling form remains valid. Our numer- 
ical data suggest that this is the case with Barkhausen 
avalanches for the range of concentrations in the interior 
region between lines (i) and {ii) (see Fig. 3). We find 
that the following scaling form is satisfied |2Q] 
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By fitting the data for L = 100, 200, and 300 to the 
expression (^ we determine the exponent a(c) and the 
fractal dimension of avalanches Ds{c) for several values 
of c (and /) in the region II. The results are also shown 
in the inset to Fig. 1 as a function of c for fixed / = 0.5. 
For each value of / a distinct family of curves is found. 
For instance, the finite-size scaling fit which is shown in 
Fig. 4 is obtained at the point c =0.1, / =1 (marked by ★ 
on the phase diagram), leading to the exponents 1+a = 
1.54 ±0.02 and Dg =1.98 ±0.02, whereas at the same 
point we find r =1.294 ± 0.003. By comparing the scal- 
ing fits for different c further c-dependence of the scaling 
function in Eq. (^ of the form 'D{x) ~ x~'^^'^^ for a: — > 
was found, with numerical values of 0(c) close to t(c). 

The scaling region is further characterized by the fol- 
lowing scaling relation |2^ 



a{c) = t(c) Ds{c) , 



(5) 



which follows from (^) by choosing sL~^= w 1 and com- 
paring to (^) . It should be stressed that numerical values 
of the exponents obtained above by the power-law and 
finite-size scaling fits satisfy (within numerical error) the 
scaling relation (j5|) for several distinct values of c in the 
region II (cf. inset to Fig. 1). 

In the critical region close to the transition line co(/) 
(line (i) in Fig. 3) the finite-size scaling form fails. 
Instead, a more general expression (^ with c-dependent 
coherence length is more appropriate, similar to the anal- 
ysis in Refs. 10,1^]. According to our results, the width 
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Sc of the critical region in the direction of the c-axis is 
rather small Sc<0.05, in contrast to the c = case along 
/-axis, where the fits to the simple scaling form (^) are 
not satisfactory in a wide range of /. This indicates that 
the scaling and nonuniversality in the region II are closely 
related to the structure of the underlying spin clusters for 

The scaling region defined via Eqs. (jj) and (^) is also 
bounded from the side of strong disorder by the line c*(/) 
(line (ii) in Fig. 3), which ends on the c-axis at the per- 
colation threshold c*(0) = 0.395. Above this line the 
power-law is practically lost and strong disorder prevents 
formation of system-size avalanches (cf. the two lower 
curves in Fig. 2). We expect that in this region the self- 
similarity persists only on a finite scale ^s(c) ^ L, with 
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similar to a subcritical cellular automaton ||2^ . Here the 
exponents r* and D* are referring to the transition point 
c*{f), and the c-dependent correlation length is expected 
to behave as ^s(c — c*) ~ (c — c*)~^= 23|. More detailed 
analysis of both critical regions as well as the precise lo- 
cation of the transition lines co{f) and c*(/) is left for a 
future study. 

By varying the driving rate Sh at fixed disorder the 
distributions D{s, Sh) appear to have different slopes. 
For larger Sh we find smaller exponents, for instance, 
at c =0.2, / = 0.5 the exponent 1+t =1.2904 ± 
0.0015 found with the driving rate Sh =0.01, becomes 
1.2575 ± 0.0017 for Sh = 0.015, and 1.239 ± 0.002 for 
Sh = 0.02. Larger Sh overdrives weaker pinning centers 
thus rendering the occurrence of larger avalanches more 
probable. This picture is in agreement with the experi- 
mental results on 1.8% Si-Fe alloy 

We have demonstrated that Barkhausen noise in di- 
luted ferromagnetic samples at low temperatures depends 
on a variety of physical parameters which can be con- 
trolled experimentally, such as strength of disorder and 
the applied driving rate. The comparison with experi- 
ments on 1.8% Si-Fe alloy 1^ where the driving rate was 
varied, is promising. Our numerical data suggest that a 
more complex behavior emerges when the degree of dis- 
order is varied. This includes a nonequilibrium phase 
transition at critical disorder co(/) and a true scaling re- 
gion with nonuniversal exponents. Further studies, both 
theoretical and experimental, will help elucidate the gen- 
eral question of conditions that result in the breakdown 
of universality. 
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FIG. 1. Fractal cluster of connected spins close to the per- 
colation threshold with areas of up spins (bright) and down 
spins (gray). Nonmagnetic defects are shown as dark points. 



FIG. 3. Schematic phase diagram showing (I) region of 
1st order phase transition, (II) scaling region, and (III) re- 
gion of strong disorder. The system ceases to percolate along 
the line (ii), while a nonequilibrium phase transition in the 
universality class of RFIM occurs along the line (i). Results 
of Fig 2 are obtained by varying concentration c along the 
dashed line. The distributions in Fig. 4 are calculated at the 
point marked by 
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FIG. 2. Double logarithmic plot of the integrated distribu- 
tion of size D{s,c) vs. size s of avalanches for fixed lattice size 
L =100 and / = 0.5 and for various values of concentration 
c= 0.0, 0.05, 0.1, 0.2, 0.3, 0.35, and 0.4 (from top to bottom), 
normalized to the first point. Inset: Exponents of distribu- 
tions of size 1 -I- r, and length 1 -f a, and fractal dimension 
Ds of avalanches plotted vs. c for fixed f—0.5 in the scaling 
region. Also plotted is 1 + Ds r. 
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FIG. 4. Double logarithmic plot of the integrated distri- 
bution of size D{s,c,L) vs. s for fixed c =0.1 and / =1.0 
and three different values of lattice sizes L =100, 200, and 
300 (top) and its finite-size scaling plot according to Eq. (^) 
(bottom). 
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